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. Theoretical tidal displacements at Goddard, MD. Crosses 
Indicate times of passes observed by Goddard laser during 
a three day period. 

. Theoretical differential tidal displacements between Goddard 
and Grand Turk Island. Crosses Indicate the times of passes 
observed simultaneously by the Goddard and Grand Turk lasers. 

. Goddard apparent station movements in height (h) , latitude (<{i) 
and longitude (X) for four 24-hour arcs and one 36-hour arc 
(identified by their 1975 day numbers). The reference arcs were 
generated using the GEM 8 geopotential model and corresponding 
force model parameters of Table 1. The standard error bars 
denote the uncertainty corresponding to the random components of 
the laser range residuals. 

. Goddard apparent station movements in height only for four 24- 
hour arcs and one 36-hour arc. The reference arcs are the same 
as those used for the results of Figure 3. The standard error 
bars denote the uncertainty corresponding to the laser range 
residuals after adjustment of station position. Errors of a 
fraction of a meter are denoted as one meter on the diagram. 

The character of the connecting lines identifies the 1975 day 
number by comparison with Figure 3. 

. Goddard apparent station movements in height only for five 24-hour 
arcs. The reference arcs were generated using the PGS 558 geopoten- 
tial model and corresponding force model parameters of Table 1. 

All standard errors were less than the 15 cm denoted on the diagram. 
The character of the connecting lines identifies the 1975 day 
numbers . 

. Comparison of tenfold theoretical vertical movements (upper panel) 
with vertical movements recovered by the dynamic method (lower 
panel). All standard errors less than 0.10 m. The character of 
the connecting lines identifies the 1975 day number. 

. Coordinate system for proof of extremum method of station-plane 
identification. Ranging is done from stations A^ and A2 in the 
station-plane to the satellite at position a. The origin is 
taken at the earth's center of mass. 
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INTRODUCTION 


The ultimate objective of this investigation is the measurement of 
the. tidal displacements of the solid earth by laser ranging to the GEOS-3 
satellite. Confirmation of earth tide theory through surface measurements 
of gravity, tilt and strain has been difficult because of the perturbing 
influences of surface discontinuities, the poor distribution of stations 
and the lack of ocean tide Information. A measurement of surface displace- 
ment by laser ranging, although not entirely Immune from such effects, 
constitutes a more direct measurement of the tidal deformations and of the 
related Love numbers h^ and 1^. The accuracy of laser ranging to satellites 
has now reached a level of between 5 and 10 cm (Vonbun, 1977) and continuing 
improvements in the dynamic models for satellites and/or the distribution of 
laser stations should ultimately lead to the detection and measurement of 
the 30 to 40 cm geometric earth tide. 

The present investigation is restricted to the analysis of NASA laser 
ranging data from three stations at Goddard Space Flight Center, Greenbelt, 
Maryland, Grand Turk Island, and Bermuda in the GEOS-3 "calibration area." 
Therefore, the necessary conditions for a purely geometric solution for 
relative station positions are not fulfilled (Escobal et al . , 1973) and 
the determinations of station positions will depend to some extent on the 
accuracy of the model for the path of the satellite. Results by Smith 
et al . (1973) for the Beacon-C satellite and a single laser station have 
shown that the fit of an orbit to a series of satellite passes rarely equals 
the quality of the laser data. Trends could be seen in the residuals showing 
departures of 2 or 3 meters from the predicted orbit. Errors in the gravity 
field, station position or other aspects of the dynamic model were suspected. 
Our approach has been to determine the effects of errors in the predicted 
orbit on the measurement of station movements and to investigate a method 
designed to minimize the effect of those errors. 


EXPECTED TIDAL DISPLACEMENTS 


In general, the measured laser statlon-to-statlon distances are 

affected both by the tidal displacements of the earth's surface and by 

the direct effect of the tidal potential on the motion of the satellite. 

% 

The Influence on the orbit of the Beacon-C Satellite, for Instance, by 
the solid-earth and ocean tides (characterized by the Love number k 2 and a 
phase lag ((>) were found by Smith, al (1973) from an analysis of the 
perturbations in the inclination of the oroit. . A subsequent fit to the 
laser range data, with this tidal effect included, (using ^ 0.245 and 

(J) ■ 3?2 derived in the previous study) showed that the in, .rred mean heights 
of the laser station from 12-hour arcs were not significantly affected. We 
are assuming in the present study, therefore, that any small errors in k 2 
and i{) in the Implementation of GEODYNE will be a second-order effect on the 
relative laser station-to-satellite distances for 24-hour arcs. 

The theoretical vertical and horizontal displacements of the laser 
stations in the calibration area due to the solid-earth tide were computed 
using subroutines NOMAN 1 and, with a small modification, NOMAN 2 (Harrison, 
1971) . The geometric earth tide in the vicinity of Goddard has a theoretical 
peak-to-peak amplitude of about 40 cm in the radial direction and less than 
5 cm in the tangential direction (Figure 1), whereas the theoretical peak- 
to-peak amplitude of differential displacements between Goddard and Grand Turk, 
for example, is 10 to 15 cm in the radial direction and less than 4 cm in the 
tangential direction (Figure 2). To a first approximation the earth tide at a 
laser tracking station can be considered constant over the few minutes of a 
satellite pass. Ideally, the time sequence of satellite passes depends only 
on the orbital period of the satellite, the rate of rotation of the earth and 
the latitude of the tracking station. In practice, problems with laser ranging 
equipment and the weather reduce the number of usable passes. 

EXPERIMENTAL RESULTS 


Two different approaches to the problem of measuring tidal movements 
of laser tracking stations were investigated. One approach, termed "the- 
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Figure 1. Theoretical tidal displacements at Goddard, MD. 
Crosses indicate times of passes observed by 
Goddard Laser during a three day period. 
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Figure 2. Theoretical differential tidal displacements between 
Goddard and Grand Turk Island. Crosses indicate the 
times of passes observed simultaneously by the Goddard 
and Grand Turk lasers. 
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dynamic method," employs 24-hour arcs as references for determining pass-to- 
pass changes in apparent station position. By this method the apparent 
station movements due to errors in the predicted satellite track as well as 
the tides have been investigated by an analysis of single-station ranging to 
GEOS-3. A second approach, termed the quasl-geometric method, attempts to 
minimize the effects of unmodelled satellite dynamics on the determination 
of tidal displacements by considering two-station simultaneous ranging to 
GEOS-3 at the precise time that the satellite passes through the plane defined 
by the two stations and the center of mass of the earth. This approach takes 
advantage of the geometrical constraints Imposed by two-station ranging and 
reduces the dependence on satellite dynamics to the prediction of only the 
distance Rq from the earth's center of mass to the satellite. 

1. Dynamic method 

Description - This method employs 24-hour arcs fitted to laser ranging 
data using the dynamic model incorporated in the NASA program GEODYNE (Martin 
and Serelis, 1975), 24-hour arcs were chosen for the Investigation because 
they were not inordinately expensive to compute, yet they are long enough to 
allow the tracking station to sample one complete tidal cycle. Each 24-hour 
arc comprises 14 or 15 revolutions of the satellite but only four or five 
passes of ranging data. 

Calcuations were carried out using force-model parameters and station 
coordinates supplied by NASA. Table 1 lists the two sets of parameters 
corresponding to the two geopotential models, GEMS (Wagner et al., 1976) and 
PGS558 (D. Smith, personal communication, 1977 and Lerch et al., 1977) used 
during this investigation. In our implementation of GEODYNE, fitting an arc 
to data from four or five satellite passes corresponds to solving for a set 
of six orbital parameters at a particular epoch and for a particular drag 
coefficient. Only the direct tidal perturbation at the GEOS-3 orbit is 
modelled in GEODYNE and not the tidal displacement of the tracking station. 

A least-squares iterative procedure was employed to compute the apparent 
position of the station with respect to the fitted arc from the laser ranging 
data taken over each single satellite pass. Two different methods were adopted: 
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TABLE 1 


SUMMARY OF FORCE MODEL PAPAMETERP A ND 

— — I— <1 I nm«'W—r »» 

STATION COORDINATES USED IN GEODYNE 


Earth Gravitational Potantlal Coafflcients 

GEM 8 

PGS 558 

(Coefficients through degree 30 and 
order 28) 


i 

1 

Gravitational constant, G (meter **3/ 
seconds **2) 

3.98601300W-14 

3.98600800Dfl4 

1 

Other perturbations s 


i 

1. Lunar gravitation applied - ratio 
of lunar mass to earth mass 

1.229997D-02 

1.229997D-02 

2. Solar gravitation applied - ratio 

of solar mass to Earth mass j 

3.329456IH05 

3.329456D4-05 

3. Gravitation applied for othac planets 

NONE 

NONE 

4. Earth tides applied - lunar and 
solar effects included 
k2 amplitude 
k2 phase angle 
k3 amplitude 

0.29 

2.50“ 

0.0 

0.29 

2.50“ 

0.0 

5, Drag applied (D65 JACCHIA 1965 static 
atmospheric density model used) 

Drag coefficient 

2.3 

Adjusted 

6. Solar radiation pressure applied 

- 1 AU solar radiation pressure 
(Newtons/meter **2) 

- Reflectivity 

- Satellite Cross Sectional Area 
(Meters *<* 2 ) 

- Satellite Mass (kilograms) 

4.500D-06 

1.500 

1.437 

3.459IH02 

4.500D-06 

1.500 

1.437 

3.459IH-02 

Goddard station position data (station 
7063) 

Coordinates - Spheroid height (meters) 
North latitude 
East longitude 

9.2900 

29“1’13V8800 

283°10'18'.'5000 

17.241 

39‘’1'13'.'3507 

283°10’19V7500 


6 


TABLE I (Cont.) 


E. Bermuda station position data 
(station 7067) 

Spheroid height (meters) 
North latitude 
East longitude 


-24.091 
32’21'13'.'7636 
295*20' 37'.'8585 


F. Grand Turk station position data 
(station 7068) 

- Spheroid height (meters) 
North latitude 
East longitude 


-19.730 
21*27 '37'.'7762 
288*52 '4'.’9584 


Earth ellipsoid - semi major axis (meters) 
Flattening 


6378155.00 

1./298.255 


6378145.00 

1./298.255 
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(X) The station was allowed to move ^ i all three coordinates by 

solving for Incremental adjustments Ax, Ay, Az from an approximate 
position X, y, z by the system of linear equations! 

2(x-x^) Ax "h Ay + 2(z~z^) Az ■ 

tl*! t»h 

where is the i observed range to the satellite, is the 1 

predicted range and x^, y^, z^ are the predicted earth fixed 

satellite coordinates, 

(2) The station was constrained to move in only the radial (height) 
direction by solving for incremental adjustments Ax, Az by the 
system of linear equations: 

2(X“X^) Ax + 2(z-z^) Az » “ P^^ 

and constraint (z/r)Ax - (x/r)Az » 0 where r is the radius to 
the station. 

These methods converge to less than 1 mm in two iterations when 
the initial position is less than 50 m in error. 

Results - During the first year of laser ranging to GEOS-3 there were a 
number of 24-hour periods during which several passes of the satellite were 
observed by the Goddard laser. Four 24-hour arcs and one 36-hour arc fitted 
to Goddard ranging data only were used as references for computing apparent 
station positions for single passes of the satellite. GEODYNE calculations 
using the GEMS geopctentlal model, appropriate force model parameters and 
coordinates (Table 1) gave, during a pass, laser ranging residuals, having a 
small random component of amplitude about 5 cm, plus a systematic component 
departing 2 to 5 m from t'< arc. We believe that the random component 
indicates the precision of the laser ranging and the systematic component 
represents the inability of the computed satellite track to fit the laser 
data. The net R.M.S. residuals for the five arcs ranged from 0.45 m to 1.96 m. 
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Apparent station movements for Goddard were computed for all satellite 
passes of the five GEODYNE arcs according to the methods described In the 
previous section. Method 1, where three-dimensional station position adjust- 
ments are allowed for each pass, gave R.M.S. variations In Goddard station 
position from pass to pass of about 8 m In height and about 11 m In latitude 
and longitude (Figure 3). In this method the systematic components of the 
laser range residuals were completely absorbed Into apparent station movements, 
leaving only a random component. The amplitude of the apparent movements 
depends on the amplitude of the systematic component of range residuals 
for the pass, on the geometry of the satellite path with respect to the laser 
station, and on the duration of tracking for each pass. In general, the 
single-pass station position Is poorly determined In a direction normal to the 
surface containing the satellite path and the station. The resulting large 
apparent movements In that direction make three-dlmanslonal station position 
determinations unreliable for the purposes of the present experiment. Method 2, 
where only changes In station height are allowed, gave significantly smaller 
apparent movements (Figure A) . Except for one short pass at a low elevation 
In arc D216/217, the four 24-hour arcs gave R.M.S. variations around 1 to 2 m 
In height. An offset In station height of about 15 m Is seen, however, for 
part of the 36-hour arc. This method does not absorb the systematic components 
of the range residuals Into station movements, yet the estimated standard 
errors on the apparent movements are, In general, less than 1 m. 

A further decrease ,ln apparent station height movement was achieved by 
Method 2 by an improvement In the GEODYNE force model. Using geopotential 
model PGS 558 and appropriate station coordinates (Table 1) , the random 
component remained the same but the overall R.M.S. residuals were reduced to 
between 0.14 m and 0.65 m for the five arcs (the 36-hour arc D207/208 was 
reduced to a 24-hour arc) . Figure 5 shows the resulting apparent station 
movements for the five arcs, plotted on an expanded vertical scale. The 
R.M.S. variation in station height for the results of Figure 5 is 0.80 m and 
the corresponding value for the theoretical tidal movements is 0.11 m. The 
detection of vertical tidal movements by this method clearly requires further 
improvements in the dynamic model for the satellite. However, the stability 
in station height is now good enough to allow the suitability of 24-hour length 
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arcs to be tested. In order to investigate the tendency of 24-hour arcs to 
absorb the real tidal movements of the laser station, tenfold-amplified 
theoretical tidal movements in height were introduced into the laser ranging 
data before computation of the GEODYNE reference arcs. Figure 6 compares the 
Induced height variations with those recovered by the method. Although move- 
ments up to 1 m are seen, they do not appear to be correlated with the input 
tides. It must be concluded that 24-hour arcs are able to absorb the geomet- 
rical tidal movements of a single tracking station and that they are therefore 
not suitable as reference arcs for measurement of the geometric tides by the 
present method. The R.M.S. amplitude of adjustments in the orbital elements 
and the drag coefficient necessary to absorb the theoretical tides were found 
to be as follows: semi-major axis, 0.09 cm; eccentricity, 0.019 x lO”^; 

inclination, 4.7 milliarcseconds; right ascension of the node, 4.9 milliarc- 
seconds; argument of perigee, 0.87 arcseconds; mean anomaly, 0.87 arcseconds; 
and drag coefficient, 0.28. These adjustments are equivalent to a movement 
of the satellite orbit in space of the order of 20 cm in both the radial and 
tangential directions. The changes in orbital elements are smaller by two to 
three orders of magnitude than the variations in orbital elements of 6-hour, 
four pass, arcs for the Beacon-C satellite due to direct tidal perturbations 
on the satellite orbit (Smith et al., 1973). Future attempts to measure the 
geometric tide by the method outlined will require the selection of longer 
reference arcs that are Incapable of absorbing the tidal movements of the 
tracking stations. 

2. Quasi-Geometric Method 


Description - This is a method for determining the radial distance from 
the earth's center of mass to each of two laser stations which are simultan- 
eously ranging on a satellite. The radial distance to the satellite is assumed 
to have been determined Independently by some other means. We use the term 
quasl-geometrlc because only the satellite radial distance must be known and 
this only for a few seconds as it crosses the plane defined by the laser 
stations and the center of the earth. This is normally the best predicted 
satellite coordinate and has the smallest rate of change; typical rates for 
GEOS-3 are 5-10 meters per second. The method is based on a concept 


Figure 6. Comparison of tenfold theoretical vertical movements (upper 
panel) with vertical movements recovered by the dynamic 
method (lower panel). All standard errors less than 0.10 m. 
The character of the connecting lines identifies the 1975 
day number. 
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uuggested by Vanicek (private communication, 1975) and Nesbo (unpublished 
manuscript j 1975) for measuring the differential earth tide (Bower, 1976) 
between two simultaneously-ranging stations. 

With reference to Figure 7 (see Appendix) it can be shown that at the 
instant that the satellite crosses the station plane (g ■ 0) there is a 
non-linear relation involving only the parameters (radius of mean station 

height), 2o) (angular separation of stations), 2h (differential radius of 
stations), Rq (radial distance to satellite), and $2 (the range distances 
from stations to the satellite) and not Involving the parameters of satellite 
position X and g. By the cosine rule the laser ranges and S2 from stations 

1 and 2 to the satellite will be given by the equations: 


« Rq^ + (R^ + h)^ - 2 Rq (R^ + h) cos g Cos (w + x) 

» R.^ + (R„ - h)^ - 2R (R - h) cos g Cos (u - x) 
^ u in U In 


( 1 ) 


Setting g * 0 and eliminating X from these equations, we obtain the following 
nonlinear relation in R^, oi, h, Rq, S^, and S2: 

sin^u) [rJ 2 (Rq^ + R^^ - h^) - - 82^} + h(S^^ - 82^)^ + 

16 (R ^ - h^)^ 

0 m (2) 

sin^o) cos^u) + cos^u) jh{2(RQ^ - R^^^^ + h^) - 8^^^ - 82^} + 

The instant that the satellite crosses the station-plane, where relation (2) 
holds, can be found either from a single pass fit to the range data by the 
program GEODYNE or by an extremum method employing the parameter h. By the 
extremum method, relation (.2), although it only holds true at g = 0, can be 
used to generate values of the parameter h for satellite position where g 0. 

It can then be shown (see Appendix) that the instant of station plane crossing 
(g = 0) is identified by the point when the computed h values reach an extremum. 
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Before proceeding further we write the nonllneat equation (?) In terns 
of the non- time varying parts of the station radii » and and the tidal 
Love number Thus: 


'll 


N, +■ 


h^tj! 


'21 


N2 + 


*'2’^2i 


(3) 


where and are the equilibrium tidal displacements of the stations in 
the radial directions, and: 


\i - <^li ^21>/2 - (Nj. + N2>/2 + h2(t^^ + x^^)l2 


‘^i “ ^*^11 ■ '^2^'^li " '^2p^^ 


(A) 


Substituting expressions (4) into (2), and neglecting tidal variations In u), 
we obtain for the 1 plane crossing a non-linear equation in Rq^, S 2£» 
Nj^, N2, h2 and as which we can represent by: 

F^(Rjj£»Sj^^,S2j^»Nj^,N2,u,h2) “ 0 (5) 


Here ^01 is the radial distance to the satellite and and $2^^ are the laser 
ranges from stations 1 and 2 to the satellite, all at the Instant the satellite 
crosses the plane. Given four or more such plane crossings, separated suffi- 
ciently in time so that the coefficients of h2 are not simply a linear combin- 
ation of the coefficients of and N2, we solve for Nj^, N2, w, and h2 by the 

Newton Raphson method. We make a first estimate of the parameters N°, N2, 
and h”, then improved values Nj^, N2, o), h2 are found by correcting the initial 
values by the amounts ANj^, AN2, Ao), and Ah2 obtained by solving the following 
system of linear equations: 

3F. 8F 3F aF. 

^ (NJ)AN^ + ^ ^ (hpih^ - 

^i%l> ®2i’ ®2’ 

i = 1 , 2 ... k where k ^ 4 
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Using the Improved values of the unknown parameters a new set of differentials 
are calculated and the procedure is repeated until convergence is achieved. 

Results - Laser ranging data from stations at Goddard, Grand Turk and 
Bermuda were examined for the presence of quasl-slmultaneous measurements to 
GEOS-3 during plane crossings in the months of July and August in 1975 and 
February in 1976. Only five usable crossings were found throughout July and 
August and none at all in February. There were many other instances of plane 
crossings but, for these, laser data was not available from both stations. 


Details regarding these five passes are listed in Table 2. The time 
shown is approximately that of the plane crossing and the columns of partial 
derivatives are with respect to the function F described earlier. The plane 
is identified by numbers referring to the stations which define the plane, 
where 1, 2 and 3 refer to Goddard, Grand Turk and Bermuda respectively. The 
columns headed dj^, 62 and d^ are the calculated equilibrium radial displace- 
ments in meters at the three stations due to the earth tide (i.e. , h2 = 1.0). 


About the times of each pass, values of RQ(t) were found from 24-hour arcs 
calculated by program GEODYNE on the basis of Goddard range measurements only 
(see Section 1, Dynamic Method). A linear equation of the form of equation (6) 
but involving the unknowr.s: 4N2* ^“23 *^^2’ obtained for 

each of the five plane crossings. The nominal values assumed for the station 
coordinates were those given in Table 1. The result of solving the five 
simultaneous equations is given in Table 3 for four cases, where AR2 and AR^ 
are the calculated differences between the true radial distances to the laser 
stations determined here and the nominal distances assumed. Similarly, Aoi^^^ 
and Ao )22 are the differences between calculated and nominal station separations 
expressed in terms of great circle distance. ARq represents the mean error 
in the predicted radial distance to the satellite for the five plane crossings. 
The result shown as Ah2 is that value of the Love number h2 which satisfies 
the five equations (i.e., h2^ = 0). 


For Solution 1 the function ^^Q(t) was assumed correct and ARq was set 
equal to zero. It is known from seismological evidence however that h2 - 0.615. 
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TABLE 3 

CALCUUTED CORRECTIONS TO Rj, Rj, wJs' ^0 

CALCULATED h2 AND STANDARD DEVIATIONS FOR FIVE PLANE 
CROSSINGS. (<*>denotes an enforced value). 


Solution 

1 

Solution 

lA 

Solution 

2 

Solution 

2A 

0 (X) 

3X 

®^o 

14.284 

0.937 

15.679 

1.090 

15.335 o(R^) 

6,179 

1.148 

-2.459 

1.069 

-2.874 

1.185 o(R^) 

1.670 

-0.088 

-0.948 • 

0.047 

-0.893 

1.107 a(R^) 

0.398 

4.656 

-0.541 

5.026 

-0.655 

6.962 cj(R ) 
0 

2.406 

3.512 

0.721 

3.702 

0.-652 

4.984 a(R^) 

1.292 

0.0* 

2.160* 

0.0* 

2.361* 








A 8@cond source of error in Rq(c) considered was that due to the 
geometric effect of the earth tide on the height of the Goddard laser. 

Since this geometric effect was not taken Into account In fitting the five 
orbits »at least part of the geometric tide would be reflected In RQ(t). The 
"trie tide for Goddard Is given by h^d^ where dj^ ly listed In Table 3. 

To dev'niine the effect of this on our results the five simultaneous equations 
were adjusted on the assumption that all of the geometric tide was reflected 
in RqCc) and the system of equations then solved as before. The results 
presented in solutions 2 and 2A are the counterparts of solutions 1 and lA 
after this adjustment is made. Note that the sum of the squares of the errors 
shown for solution 2A Is slightly larger than for lA but h 2 is closer to the 
theoretical value. Presumably the correct solution would fall somewhere 
between solutions lA and 2A. 

Estimates of the effect of random errors in RqCI) on the solutions are 
presented in column 6 in terras of the standard deviation of RqCc) about true 
values. If we suppose either that the nominal station coordinates adopted 
for this analysis are correct within 1 - 2 m or that the Love number is known 
to be 0.6 then the results suggest that RQ(t) is systematically less than the 
true values by about 2.0 m with a much smaller random error. 


CONCLUSIONS 


The 5 cm precision of laser ranging measurement i® certainly adequate 
for the observation of the 40 cm geometric earth tide, or even the 15 cm 
differential tide between two stations which can observe a satellite simul- 
taneously. However, we cannot yet predict an orbit based on one tracking 
station and 24 hours of data which is stable enough to be used as a platform 
to observe the total tidal displacements. 


The present dynamic model, employing gravitation field model PGS 558, 
fits 24-hours of laser data from a single station leaving systematic residuals 
during a pass of up to 1 meter, and resulting in apparent station movements in 
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height o£ comparable amplitude, thus hiding the tidal variation. But even in 
the absence of Imperfections in the dynamic model, 24-hour arcs would tend to 
absorb the tidal movements of a single tracking station leaving the pass-to- 
pass apparent station heights unchanged. Longer arcs or arcs using data from 
more than one laser should be less likely to absorb the geometric tides, but 
they would be expected to fit the laser data less well. 

The quasi-geometric method is influenced significantly less than the 
dynamic method by errors in the predicted satellite position because this 
method only requires a knowledge of the radius to the satellite and it is 
8en?ietve principally to the differential tidal displacements between laser 
stations. Due to the stringent conditions that must be met for a usable 
plane-crossing to occur, however, there has been difficulty in finding a 
sufficient number of plane crossings for a rigorous statistical test of this 
method as it is presently implemented. But, for five passes over the cali- 
bration area that satisfy the criteria, a good approximation to the theoretical 
Love number h^ is obtained when a systematic bias of 2.16 meters is allowed 
in the radial distance to the satellite. This bias is justified independently 
by the assumption that the nominal station coordinates are correct within 
1 - 2 ra. The value of h 2 appears to be reasonably insensitive to changes in 
the predicted radial distance to the satellite due to absorption of the tidal 
movements of Goddard by the GEODYNE reference arcs. 
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APPENDIX: PROOF OF THE EXTREMUM METHOD OF STATION-PLANE IDENTIFICATION 


Let us consider a reference system of spherical coordinates with origin 
at the center of the earth 0 and axis OZ perpendicular to the plane Q^Aj^OA2p2 
passing through the center of the earth. The radial directions to the 
ground stations A^ and A2 are extended to points and Q2 such that OQ^ > 

OQ2 ■ Rq, the radial distance of the satellite at any Instant from the center 
of the earth. T' axis of meridional reference OX In this plane bisects the 
angle Qj^UQ2* Then, with g as the perpendicular arc from the position 0, of 
the satellite to this plane and x as the arc from the foot of this perpen- 
dicular to the bisector OX, the coordinates of tne satellite can be denoted 
as (Rq, g, x). Also, A^ and A2 can be represented In the same system of 
reference by the coordinates (R^^ + h, 0 , -w) and (R^ - h, 0 , m) respectively, 
where R^ 1 ^ the mean radius to the ground stations, 2 h Is their elevation 
difference and 2 o) Is the angle they subtend at the center of the earth. Then 
laser ranges and $2 which are linear distances from the satellite 0 to 
ground stations A^ and A2 respectively will be given by the equations 


= Rq^ + (R^ + h)^ - 2 Rq(R^ + h) cosg cos(oj + x) 
® 2 ^ + Rq^ + (R^ - h)^ - 2 Rq(R^ - h) cosg cos(o) - x) 


(Al) 


These ate the basic equations which are used to evaluate h from the 
observations and S2. In developing our method of solution for h, we 
assume that w and R^^ are constants which are known before-hand. Range 
observations and S2 as well ns the radial distance to the satellite R^ 
which vary with time are assumed to be available at discrete instants of time. 


The problem is not solvable in its present form since, corresponding 
to n given sets of (S^^, S2, Rq) values, we have n pairs of equations of the 
type (Al) involving 2 N + 1 unknowns viz., n difference g's, n different x’s 
(as g and x vary with time) and a constant h. Thus, the number of unknowns 
being more than the number of equations by one, no unique solution for h 
will be possible until and unless an additional condition for the problem 
is made available. 
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To obtain this, we consider a second pair o£ equations 


where H, unlike h, and X, different from x, vary with time. 


(A2) 


Eliminating X between the two equations in (A2), we can write 

16R ^(R„^ " - sin^a)lR^{2(RQ^ ^ ‘ " ®2^^ ■*■ 

0 ' n 

H(Si^ ~ + cos^ui[H{2(Rq^ " 

(A3) 

from which H can be obtained when other quantities are known. We have 
developed a subroutine which computes H iteratively, starting from the initial 
value of H = 0. 

If we assume the x-eliminate between the equations in (Al) can be 
formally written as 


h = F(Sj^, S2 » Rq, cosg) (A4) 

then, the similar equation for H will be 

H = F(S^, S^, Rq, 1) (A5) 

which is another form of (A3) . 

Differentiating (A4) with respect to time t and remembering that h is 
independent of t, we have 


0 . ^ (Si.Sj.Rg.coag) • ^(Sj.Sj.Ro.cosg) ' ^ 

vS^ X ^^2 " 


3F 

9R« 


(Sj.Sj.Ro.cosg) • ^ ^RCSj.Sj.Rj.cosg) j ' 


(A6) 


0 dcosg 


Substituting In (A6) t “ t^ corresponding to g « 0, we have 
0 . 3F(S^rS2,Ro,1) • ^^(S,.S2,Ro.1) ' ^(S^rS^.Rq.I) ‘ 


(A7) 


dS 


1 9S 


2 ’ 3R 


Q 


which, when compared with the equation obtainable from differentiation of 
(AS) with respect to t, yields: 


H ® 0 when g = 0. 


(A8) 


Consequently, from comparison of (Al) and (A2), we find that when 
H = 0 corresponding to g * 0, H = h. 


Thus, the equation (A8) Is the additional relation th&t Is needed to 
obtain h uniquely. 


In practical computation, values of H are computed Iteratively from each 
set of (Sj^jS^.Rq) values, using the subroutine based on (A3). The plot of 
these values of H against time would show a smooth curve with an extremum 
(l.e. , H = 0) occurring at the Instant when the satellite crosses the vertical 
plane through the ground stations (l.e. g = 0) . For precise computation of 
this instant, four consecutive values of H are selected such that 
and H(tj^^ 2 ) either both greater than or both less than H(tj^) and 
A third-degree polynomial in time is then fitted to these values of H to obtain 
the extremum value for H and this is the value of h we require. This part of 
the computation is implemented by a second subroutine. 
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